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ABSTRACT 

The role of a resonant z^e-^^s osciUation is discussed in the event of a supernova explosion. 
It is concluded that a significant Ue-^s mixing may hinder the ability of the supernova to 
explode. It may also cool the proto-neutron star too quickly with respect to the observed 
cooling time of several seconds. The constraints on the Ve-T^s mixing parameters based on 
the above arguments are calculated. 
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1. Introduction 



Oscillations between active neutrinos v^, and sterile neutrinos Vs have been discussed in 
various astrophysical and cosmological environments. Such examples include the oscillation 
between i/g and Vs as solutions to the solar neutrino problem [1] , and oscillations between 
Vq, and Vg in Big Bang Nucleosynthesis which may affect the primordial ^He abundance 
[2-6]. Constraints on neutrino mixing parameters have been obtained by requiring the 
predicted ^He abundance under these neutrino mixings not to exceed the observed value 
[3-5]. Sterile neutrinos are also suggested as dark matter candidates if they have proper 
mixings with active neutrinos [6]. In this paper we attempt to consider Ve-^s mixing in 
another astrophysical environment-the Type II supernova (hereafter we omit "type II" for 
simplicity) core. 

Since there isn't a fully convincing numerical supernova model at present, we try to 
construct our arguments as model-insensitively as possible. Supernova explosions occur 
when massive stars (> SM©) reach the end of their evolution [7]. For such massive stars, 
an iron core is developed at the center of the progenitor star, surrounded sequentially by 
shells of silicon, oxygen, nitrogen, carbon, helium, and finally an envelope of hydrogen. 
Nuclear reactions occur at the interfaces of these shells and between the shell and the 
core. The initial iron core of the progenitor star has a density of about 10^ g/cm^, and 
is primarily supported by the pressure of the degenerate electron gas. As the nuclear 
burning of the silicon shell continues, more iron is deposited into the core, until it reaches 
the Chandrasekhar mass IAMq)^ when the core cannot be sufficiently supported by 
the electron degeneracy pressure and begins to collapse. In the meantime, due to the high 
density and the high temperature of the core, the deleptonization process e~ -\-p ^ n-\-Ve 
occurs rapidly. When the density reaches lO"'^-'^ g/cm^, the mean free path of neutrinos I 
(fa 10^/pi4 cm for neutrinos of ~ 20 MeV, where pi4 is the core density divided by IQ-*^^ 
g/ cm^) is shorter than the radius of the core -Rcore ~ 10— 100km. When the density is higher 
than 10^^ g/cm^, the diffusion timescale of neutrinos becomes longer than the timescale 
of the collapse of the core. Neutrinos are thereafter trapped inside the core instead of 
escaping freely. The collapse of the core finally stalls and bounces back when the center of 
the core reaches approximately nuclear density fewxlO^'^ g/cm^, at which the pressure of 
the degenerate neutron gas stops the gravitational collapse. The bounce creates a shock 
wave in the midst of the core that propagates outwards through the infalling matter. It 
is generally believed that if the core is small (< IAMq)^ the shock wave will have enough 
energy to blow the matter outside the core away and lead to an explosion event [8] . If the 
core is sufficiently massive, however, the shock wave quickly stalls within the core by losing 
energy in dissociating large quantities of infalling nuclei into free nucleons. It is currently 



2 



believed that in this case neutrinos deposit their energies into the shock several tenths of 
a second after the bounce and reheat the shock wave to eventually blow the outer layer 
away [9]. 

It is estimated that about 99% of the energy is emitted in the form of neutrino bursts. 
Therefore, the role of neutrinos in the supernova event is essential [9]. At the end of the 
collapse, electron neutrinos are trapped instead of streaming out of the core. The i^^, v^-, 
and anti-neutrinos are also produced thermally inside the hot core with an average energy 
of several tens of MeV. This leads to a neutrino sphere with an average energy of ~ 20 
MeV trapped inside a dense core of nucleons and heavy nuclei. The number density of 
Pe is highly suppressed at the collapsing stage of the supernova due to the existence of 
the degenerate v^. gas. Its number density builds up only after significant amount of Vf, 
escapes from the core after deleptonization. It was shown by Wilson et aJ. that the v and 
V annihilating at the neutrino sphere can deposit their energies in the shock and lift the 
shock front outwards, thus leading to an explosion [10]. Diffusive neutrinos further cool 
down the core (a proto-neutron star) to a neutron star after the explosion, on a timescale 
of about 10 seconds. This diffusive Pg Aux has been observed by Kamioka [11] and 1MB 
[12] in conjunction with SN1987A, with a typical energy of 15-20 MeV and a timescale 
of about 10 seconds. Therefore, the current picture of a supernova event has received 
significant observational support. 



The above picture assumes the standard electroweak theory. If neutrino mixing is 
present, the picture may change substantially. Many literatures have been focusing on the 
constraints on properties of neutrinos [13]. We will focus on the impact of a Ue-^s mixing 
on the supernova event, which hasn't been explored in full detail before [14]. In particular, 
instead of discussing the general cases of the mixing, we will discuss only the case of a 
resonant oscillation between and z/^, which are relatively less model-dependent. 

In the flavor basis the effective Hamiltonian of the (active)-z^s (sterile) mixing system 
in a medium is given by [15] 



where E is the neutrino energy and = m\ — m\ is the difference of the squares of the 
two vacuum mass eigenvalues mi and m2. We assume > when the sterile neutrino 
is heavier. Furthermore, 6 is the vaccum mixing angle and V^^ is the effective potential of 
the active neutrino v^, in the medium. In a stellar environment [16], 



2. Neutrino Mixing in the Supernova Core 
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= %/2GF(-0.5iV^ + iV,^ + 2N,^ + iV,J (2) 

where N^, and AT;^^ are the number densities of electrons, neutrons, and neutrinos minus 
the number densities of their anti-particles, respectively. V,^^ is similar to V^^^ except that 
and are interchanged in the equation. The effective potential of a sterile neutrino is 
always zero. Based on CP symmetry, Vj^^ = —Vj^^. In the special situation of a supernova 
core, N,^^ and N,^^ are negligible, since and v^- are only produced in vV pairs via neutral 
current processes. 

A resonance between the two mixing neutrinos occurs when the diagonal elements of 
the effective Hamiltonian eq. (1) become equal, i.e., when the neutrino energy equals the 
resonance energy Ej. given by 

^ dm? cos 29 

Er = . (3) 

Inside a supernova core, where spherical symmetry is a reliable approximation, the region 
where such a resonance occurs for neutrinos with an energy Er is a spherical shell with 
radius Rres{Er)- It can be shown that if the sterile neutrinos resonantly converted from 
the active neutrinos escape freely from the core (we will discuss the condition for this 
later), the usual continuity equation relating the local rate of change of the active neutrino 
number density n^^ to the divergence of the corresponding local fliix (p^^ is modified to 



n,, = -V • 



iout j^out 



p.,- (4) 



Here, (p°^^ and n°^* are the partial flux and number density which only account for z/q,'s 
going outward with respect to the surface of constant resonance energy E^ (outward radial 
direction in our case). Furthermore, is the probability for an active neutrino going 
through the resonance to be converted to the sterile neutrino Us- We have neglected terms 
corresponding to z/g z/q, convertions and z/q, — > Ug convertions of inward going z/q,, since 
we are only interested in the global z/q, loss by integrating eq. (4) over the whole SN core 
to which the sum of the neglected terms does not contribute. 

In general, Pq depends on the rate of decoherence of the neutrino flavor density matrix 
caused by non-forward scattering off the background [14] as well as on the adiabaticity of 
the resonance. In the limiting case where the width of the resonance region is much shorter 
than the mean free path of neutrinos, which is the case we will do our calculations. 



-TrJm^ sin^ 26 



E 



-1 



Pa{Er) 1 - exp VEr + — . (5) 

4 cos 20 c 

4 



This is correct to the first order in the change of Er in the resonance region [15]. Pq ~ 1 
corresponds to a highly adiabatic resonance whereas Pq, <^ 1 represents a highly non- 
adiabatic resonance. On the other hand, when the width of the resonance region is much 
larger than the neutrino mean free path, ~ 0.5 [14], which means the resonance is 
always "adiabatic". Therefore, it can be seen from our following calculations that our 
constraints that exclude sufficiently adiabatic resonant mixings will equally apply to the 
case with large resonance width, as long as the Vs stream freely outside of the resonance 
region. 

To get an order of magnitude estimate of the neutrino loss rate, we can further simplify 
eq. (4) by setting VEr ■ (d0°^Vc/E) (K) ~ idEr/dr)ncf(Er)E^/E^ and {dn^'^y dE){Er) ^ 
nf{Er)E^/E (we drop the subscript Ua from now on). Here, c is the speed of light, f{Er) 
is the occupation number at a state with Er, and E is the average neutrino energy which 
for the partially degenerate Ue is of the order of its chemical potential /i and for other 
neutrino species (and all antineutrinos) is 3 times their temperatures in the supernova 
core. We can then integrate the second term on the r.h.s. of eq. (4) over the SN core to 
get the additional global Uq, loss rate iVres and the additional global energy loss rate Q^es 
due to the resonant conversion, 



Nr. 



Qr 



I 



Rc 



Rc 



dE. 



E. 



Airr^n—^ ic + Vr) f{Er)^P{Er)dr, 
dr jt; 



' A7rr''n^{c + Vr)fiEr 
dr 



Ej, 



P{Er)dr, 



(6) 



where Vr = Er/\VEr\, the speed of the resonance shell for neutrinos with energy E^. 
The integrals on the r.h.s. are dominated by integrands at E^ ^ E if the resonance for 
neutrinos with energy E ever occurs inside the core. When E^ ~ E, we always expect 
\E\ « c ■ dE /dr since otherwise the neutrinos will build up inside or dissipate from the 
core in Rcore/c ~ 10~^s, which contradicts our knowledge of supernovae. Therefore Vr and 
E are negleted in the following discussion. If we denote A^tot and Qtot as the total number 
and the total energy of i>a inside the core, we can approximate eq. (6) by 



iVres ~ 



T-J ~ tot 7 ^ i T-J 

-'•'Core -'I'core 

One crucial quantity in the problem is the scale height of E^, \dliiEr/dr\~^ . A rough 
estimate for it is 

d]-nP,^ -1 d]r,N^^ -1 HAun -1 N^^ 

(8) 



Qr 



cP{E)\ 



Er=0 



Q 



tot- 
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dlnEr 


-1 


cJlniVeff 


-1 


dlnn 


-1 


iVeff 


dr 




dr 




dr 




n 



where \dliin/dr\ ^ is the scale height of the neutrino density n at the resonance region, 

and A^eff = Vi^^/^^F- Without going into detailed models, 

dlnn 



dr 



> I. 



(9) 



In the following sections we apply these results to estimate the loss rates in the col- 
lapsing and the cooling phases of supernovae. 

3. Constraints from the Collapsing Phase 

During the collapsing phase of supernovae, Vg. will encounter a resonance when the 
Ve-^s mixing satisfies 

5m^ cos 26 = 2V2GFErNe{{ = 2V2GFEr{Ne - 0.5Nn + 2N^J. (10) 

Taking E ~ lOMeV, and the density of the core p ~ 10^^ g/cm^ to several times 10^^ g/cm^, 
we find 10^ eV^ < 5m? ^10^ eV^. At this stage the VeS are already trapped in the SN 
core. The relevant time scale here is the free fall time scale of about 10~^s~ i?core/(0.1c). 
Then if Nres > (0.1c/i?core)^tot) J^e wiU never get trapped after the central density of the 
core reaches and exceeds the resonant value. The process is illustrated in Figure 1. It 
has been shown that no explosion occurs when the trapped lepton number is less than 0.9 
times the lepton number before the collapse [17]. Therefore we require 

P(E)\^^Q<0.1 for lO^eV^ < < loW (11) 

so that no serious leakage of will occur. Because of eq. (5) this corresponds to 

-(^)"'^0.1. (12) 



5m'^ sin^ 29 


dlnEj. 


E 


dr 



Since iVeff- ~ n, \d\nEr/ dr{E)\-^ > I ~ lO^cm from eqs. (8) and (9). Taking E ~ lOMeV, 
we conclude that Ue-i's mixings with 

sin^ 2^ > lO-^eV^ and lO^eV^ < ^m^ < lO^eV^ (13) 

are ruled out so that the supernova core can trap enough lepton number to explode. 

4. Constraints from the Cooling Phase 

For < 10^ eV^, interesting constraints on the Ug-'^s mixing come from the con- 
sideration of the resonance that occurs after the bounce. After the bounce, N^ff = 
Ne — O.bNn + 2Njj^ drops below zero very quickly at the outer part of the core [17,18]. (In 
general = Np ^ before the collapse, where Np is the number density of protons. 
After about 1/3 of protons are deleptonized into neutrons, A'^gff will be negative.) However, 
at the center of the core, there is a significant build-up of z/g, N^f^ remains positive. Figure 
2 shows roughly the relation between A'^efj and the position inside the core at about 0.5 
seconds after the bounce, as inferred from figure lOf of ref. 18 (similar relations exist for 
prompt-explosion models [17]). Thus, there exists a central part of the core with a positive 
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Vi,^ and an outer part of the core with a negative V,y^. Furthermore, outside the core where 
no significant deleptonization has taken place, Vi,^ is still positive. Therefore two resonance 
regions may form, (for both i/g-^s mixings and Pg-^s mixings), one inside the core and one 
outside the core, with 

iVeff ~ 0, (14) 

if 5m^ is much smaller than GpENQfi and can be neglected. If we take ~ 10 MeV at 
the edge of the neutrino sphere, we can see from figure 2 that when Sm^ < 10^ eV^, both 
Ve and Pe can be resonantly converted into sterile neutrinos inside the core. 

In this section we consider the resonance region inside the core. In this case, the scale 
height of Er is much smaller than in the previous section due to eqs. (8) and (14). By 
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using eq. (10) and n < E /Stt^, we find for Er ~ E, 



dlnEr 



dr 



-1 



n 



{E)>r-^ > I ^ (15) 



^/2GyE 



The concern is that if the resonance occurs with a sufficient adiabaticity, the intensive 
production of the freely streaming sterile neutrinos will cool down the core too quickly 
compared with the observed cooling time of about 10 seconds. This translates into 

cP\ _ — _ 

"(^)>lsec-\ or equivalently ^^{E) >1Q-'^ (16) 



\Er.=0 



for -Rcore ~ 10 to 100km. By using an approximation 

— 5rn? siv? 29 



P\E.=oiE)- ^^^^Q |VK|-^ forP«l, (17) 



we conclude after combining eqs. (15), (16) and (17) that if 

^ cos 2^ ) y^G^E' 

the proto-neutron star will cool down much faster than the observed timescale by losing 
neutrinos. Taking E 50 MeV, / ~ 10"^cm, eq. (18) excludes Ug-Ug mixings with 

(5m^sin2^ > 30eV^ and < lO^eV^. (19) 

5. Possible Tightening of the Bound by Reconversion 

Since the effective potential V^^ goes from negative at the edge of the core to positive 
outside the core, at least the Ue-i^s will encounter a second resonance if < 10^ eV^. 
The density outside the core is quite low compared with that of the core. As a result, 
this outer resonance could be highly adiabatic since \dliiEr/dr\~^ could be much larger. 



Then all the Vg converted from deep in the core will be completely reconverted into 
which then freely escape from the supernova because the outer resonance is located above 
the neutrino sphere. But these neutrinos have energies >50 MeV as they are originated 
from Pe's deep inside the neutrino sphere. Since the cross section of neutrinos in the v-e 
scattering experiments (Kamiokande and 1MB) scales roughly with where E is the 
energy of the neutrinos, the fact that we didn't see any V^. with > 50 MeV in these 
experiments will put a tighter bound on Ve-^s mixings than that of the previous section. 
The average energy of Pe observed by Kamiokande and 1MB is ~ 20 MeV [13]. Therefore, 
the above consideration could extend eq. (16) by two order of magnitude, because even if 
the energy loss rate due to the resonant conversion is only one tenth of the standard loss 
rate, ~ 5 of the 20 neutrinos seen by 1MB and Kamiokande should have had an energy 
of >50 MeV. More generally, if in the future a neutrino detector detects M supernova 
neutrinos with none of them having an energy > 50 MeV, the excluded parameter region 
would extend to 

5m'^sui2e > (^^^ eV^ and Sm^ < lO^eV^. (20) 

With the current data, jV = 20, the excluded region is 

(5m^sin2^ > leV^ and < lO^eV^. (21) 

One has to keep in mind, however, that eqs. (20) and (21) are only true if the outer 
resonance is highly adiabatic, which for eq. (21) requires the density scale height at the 
outer resonance to be larger than ~ (10/sin2^)m, which seems fairly reasonable when 
compared with the core dimensions. 

6. Discussion and Summary 

As we mentioned in the previous section, the above calculations are valid only if the 
sterile neutrinos have a mean free path larger than the radius of the core after passing 
the resonance. Outside the resonance region, the mean free path of sterile neutrinos is 
rouK hly // sin^ 9^ [3-5,14], where 9 TO is a typical medium mixing angle. In the case of the 
constraints from the collapsing phase discussed in section 3 this is of the same order as the 
vacuum mixing angle 9 whereas in the case of section 4 due to medium effects this is even 
suppressed compared to 9. Our calculations are thus valid as long as 

sin^6' < Z/i?core ~ 10"^ (22) 

Figure 3 shows the excluded region on the z/g-z/g mixing parameter space from our 
consideration of a resonant Ue-Ug mixing in both the pre-bounce and the post-bounce 
supernova core, eqs. (13) and (19), as well as the possible tighter bound eq. (21), and 
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eq. (22) under which the above bounds are vahd. It also shows the parameter space 
excluded from Big Bang Nucleosynthesis, which gives a stricter constraints for Sm^ < 10^ 
eV^ [3-5]. Tha hatched region in Fig. 3 shows the required mixing parameters for the 
sterile neutrino to contribute between 30% and 100% of the critical density of the universe 
[5,6] (assuming a Hubble constant between 50 km/sec/Mpc and 100 km/sec/Mpc). A 
significant part of the hatched region (at which the sterile neutrino serves as warm dark 
matter) is excluded by our constraints from the supernova. 

In summary, we have considered the effect of a resonant I'e-J^s mixing on a type 
II supernova. By requiring that a supernova retains enough leptons to explode and is 
not cooled down within a fraction of a second, we obtain constraints on the Vg-J^s mix- 
ing parameter space. The excluded region includes the region required for the sterile 
neutrino to be the warm dark matter through Vg-i^s oscillations in the early universe. 
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Figure Caption: 

Figure 1. An illustration of the configuration of inside the core during the core collapse 
with (the dashed line) and without (the solid line) the Ue-i'g mixing that satisfies eq. (13). 

Figure 2. An illustration of A^gff inside the core at the bounce and several tenths of a 
second after the bounce. Rns is the radius of the neutrino sphere. 

Figure 3. On the Vg-J^s mixing paramter space, the area enclosed by the solid line is 
excluded by our supernova consideration. The area enclosed by the long-dashed line is 
the possible extended bound due to the resonance outside the core. The area enclosed by 
short-dashed line is excluded by Big Bang Nucleosynthesis. The hatched region shows the 
parameters for a sterile neutrino dark matter candidate (generated a Vg-^^s mixing in the 
early universe) to contribute 30% to 100% the critical density of the universe today (the 
Hubble constant is taken to be between 50-lOOkm/sec/Mpc). 
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